function y=Kaifang(a,eps,x0)
%a是被开方数
%eps表示精度指标
%x0表示初值
%y是a的开方
x(1)=x0;
x(2)=(x(1)+a/x(1))/2;
k=2;
while abs(x(k)-x(k-1))>eps
    x(k+1)=(x(k)+a/x(k))/2;
    k=k+1;
end
y=x';